function [CL, CD, alpha, alphadeg] = getCLCD(V, Vc, N, deg2rad)
    CL0 = -0.041065;
    CL1 = 0.016292;
    CL2 = 0.0002602;
    CD0 = 0.060505;
    CD1 = -0.03026;
    CD2 = 0.86495;
    alphadeg = nan(1, N+1);
    alpha = nan(1, N+1);
    CL = nan(1, N+1);
    CD = nan(1, N+1);
    for j = 1:N+1
        if V(j)*Vc > 4570
            alphadeg(j) = 40;
        else
            alphadeg(j) = 40 - 0.20705*(V(j)*Vc - 4570)^2/340^2;
        end
        alpha(j) = alphadeg(j)*deg2rad;
        CL(j) = CL0 + CL1*alphadeg(j) + CL2*alphadeg(j)^2;
        CD(j) = CD0 + CD1*CL(j) + CD2*CL(j)^2;
    end
end